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We derive exact analytical solutions describing multi-soliton complexes and their interactions on 
top of a multi-component background in media with self-focusing or self-defocusing Kerr-like non- 
linearities. These results are illustrated by numerical examples which demonstrate soliton collisions 
and field decomposition between localized and radiation modes. 
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I. INTRODUCTION 

The multi-soliton complex (MSG) is a new concept of 
modern nonlinear science If we think of solitons as 
elementary particles, then we can imagine multi-soliton 
complexes (MSCs) as combinations of elementary parti- 
cles such as atoms. In most cases an atom is an isolated 
particle, but it can demonstrate more complicated behav- 
ior. Mathematically speaking, a multi-soliton complex is 
a self-localized state which is a nonlinear superposition of 
several fundamental solitons. As it is a composite struc- 
ture, we can expect it to behave in a more complicated 
way than a single fundamental soliton. At the same time, 
an MSG can behave as a single particle unless there is a 
substantial external perturbation which causes it to split 
into its primary constituents. 

There are various phenomena in nonlinear physics 
which can be treated as MSGs, some of which are listed 
in Ref. 0|. Here, we are mainly interested in optical ap- 
plications. In fiber optics, an MSG is a single short pulse 
formed by a nonlinear superposition of fundamental soli- 
tons trapped in a common potential well, so that all of 
them have the same velocity . The nonlinear superpo- 
sition can be either coherent or incoherent, in the sense 
that the phases of the separate solitons in the bunch can 
be either related or independent. The solitons in the 
group may be bound together, or at least may stay close 
to each other, simply because they have the same ini- 
tial speed. The latter happens in the case of integrable 
systems. 

Incoherent solitons in space 10,^ have recently at- 
tracted considerable attention [bT, especially after the 
first experimental demonstration of their existence 
Incoherent self-trapping in a biased photorefractive crys- 
tal is usually well-described by a set of M coupled nonlin- 
ear Schrodinger equations (NLSEs) with saturable non- 
linearities Moreover, in photorefractive media with 
a drift mechanism of nonlinear response, the Kerr-like 
nonlinearity is a good model 0. In the latter case, the 
equations become integrable, and this allows us to obtain 
solutions analytically. Moreover, the solution can be con- 
sidered as a nonlinear superposition of a certain number 
of fundamental solitons and radiative waves. This fact 
allows a simple qualitative approach to the problem, in 
addition to having an exact solution. 

The number of solitons in a complex depends on the 
nature of physical system, as well as on the initial con- 



ditions. Mathematically speaking, we will talk about N 
fundamental solitons all having the same speed and mov- 
ing (or resting) as a single complex. The value of N can 
be 1 (fundamental soliton) but can take large values up 
to infinity. 

Elementary particles in nature never exist in isolation, 
but are submerged in a fluctuating radiation field. Simi- 
larly, MSGs can exist on top of a background plane wave 
and, more generally, on top of an arbitrary number of 
plane waves propagating in different directions. In this 
sense, we will talk about a sea of radiation modes, which 
can be stable with respect to modulational perturbations 
in both self- focusing and defocusing media Plane 
waves can also serve as a basis for the existence of dark 
solitons. An incoherent soliton is an example of an MSG, 
and stable propagation of both dark and bright spatial 
optical solitons on an incoherent background has been 
observed in recent experiments P,p^. 

In this work we present a class of exact solutions to 
the coupled set of M equations which describe symmetric 
and asymmetric MSGs on a background which is formed 
from a nonlinear superposition of radiation modes. This 
differs from the earlier results presented in Refs. [pd 12 
where the background consisted of a single mode, and 
Refs. [pp^, where only some stationary solutions with 
sech^-type intensity profiles were considered. Our new 
exact solution of MSGs in a sea of radiation modes ac- 
counts for an arbitrary number of both soliton and ra- 
diative components. 



II. MODEL EQUATIONS AND THEIR 
PROPERTIES 

It can be shown that, for photorefractive media with 
a drift mechanism of nonlinear response, a good approxi- 
mate model describing the propagation of M self-trapped 
mutually incoherent wave packets in a planar waveguide 
is the set of NLSEs for a Kerr-type nonlinearity 
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where V'm denotes the m-th component of the beam, x 
is the transverse coordinate, z is the coordinate along 
the direction of propagation, a is the coefficient repre- 
senting the strength of nonlinearity, and dn is the change 
in refractive index profile created by all the incoherent 
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components in the Ught beam. Because the response 
time of the nonUnearity is assumed to be long compared 
with temporal variations of the relative phases of all the 
components, the medium responds to the average light 
intensity, and this is just a simple sum of modal inten- 
sities [Q. Then, after defining a new set of functions 
Um{x,z) = ^/\a\tpm{x, z), wc rcducc Eq. to the fol- 
lowing set of normalized equations: 
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Here 1 < to < Af, M is the number of components, 
while s = in a self- focusing medium and s = —1 in a 
defocusing material. 

Let us briefly outline some general properties of 
Eqs. (^). The most remarkable fact is that these cou- 
pled equations are completely inteqrahl e by means of the 
inverse scattering technique (1ST) ||l4|-[lq|. Specifically, 
any solution can be represented as a nonlinear superpo- 
sition of solitary waves and radiation modes which cor- 
respond to the discrete and continuous parts of the 1ST 
spectrum, respectively. 

Every fundamental soliton (labeled j) is character- 
ized by several eigenvalues: (i) a complex wavenumber 



i/ij, (ii) a shift in the coordinate plane 
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and (iii) a polarization vector p^-i-* in the functional space. 



normalized to unity as ^ 



M 



U) 



1. The simplest 



bright single-soliton solution in a self-focusing medium 
(s = -t-1) can be written as: 



ix,z) =Pm rj sech(/3j) e''^^ 
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where (3j = r^ {xj — fijZj), = fijXj + {rj — fi'j)zj/2, and 
i^Xj ^ Zj ) — {^X X j , Z Zj ) are the shifted coordinates. We 
see that the peak soliton intensity and its inverse width 
are determined by the real part of the wavenumber r-,-, 
while the imaginary part /ij characterizes the soliton ve- 
locity in the transverse direction. Moreover, each funda- 
mental soliton can be "spread out" into several incoher- 
ent components, as defined by the polarization vector. 

The solution for a single radiation mode, in the form 
of a plane wave, can also be characterized by a similar 
set of parameters. 
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where CFj — fixj + (2sr| — /i^)zj/2. Such a plane wave 
exists for either sign of nonlinearity, s = ±1, and it is 
stable in a self-defocusing medium. Moreover, an inco- 
herent superposition of a large number of plane waves 
can be stable even in the self-focusing medium, as was 
shown in 1^. In the presence of solitons, the plane waves 
are distorted, but due to the integrability of the original 
equations, the corresponding solutions can be obtained 
in an explicit form, as we demonstrate in the following 
section. 



III. EXPLICIT SOLUTIONS FOR 
MULTI-SOLITON COMPLEXES ON A 
BACKGROUND 



It has been demonstrated in earlier studies |12 17 
that a stationary MSC can only be formed by in- 
coherently coupled fundamental solitons and radiation 
waves which have orthogonal polarization vectors, i.e. 

X]m=iP™V™^ = '^i"- We restrict our investigation to 
this important case. Then, the mathematical description 
can be simplified if we use the rotational symmetry of 
the functional space of the original Eqs. (^). Indeed, it 
is sufficient to find solutions Uj where all the fundamen- 
tal nonlinear cigenmodcs belong to different components, 
Pni = Smj, and then the full family of solutions can be 
determined using the following transformation: 
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where the matrix Rmj defines a rotation in the M- 
dimensional space (characterized by M—l angles), 
which preserves the MSC intensity profile J2m l^mP = 

As follows from the 1ST, it is possible to construct 
the full solution by adding the radiation modes to bright 
MSCs. On the other hand, solutions for bright MSCs, 
having a zero far-field asymptotic in the self-focusing 
medium (s = -fl) can be found by solving a set of linear 
equations |18 1^: 
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where Ms is the number of fundamental solitons, and 
Cj = Xj 6xp {(3j + i"fj ) . Here the coefficients Xj define the 
relative coordinate centers for the fundamental solitons. 
By choosing their values in a special way, the solutions 
of Eqs. (^ can be presented in a simple, symmetric form 
(see Refs. [^^^ for details). 

In order to reveal the basic properties of radiation 
modes under the presence of an MSC, we perform a linear 
analysis, assuming that the wave amplitudes are vanish- 
ingly small and that they do not contribute to the inten- 
sity profile. We note that Eqs. (^) define an extension 
of the functional basis introduced by Kay and Moses |2^] 
to construct reflectionless potentials. The difference is 
that in Ref. |21| the kj are real, whereas we performed a 
similar analysis and found that the results can be gener- 
alized to the case of complex wavenumbers. Specifically, 
we derived solutions for scattered plane waves, or radia- 
tion modes belonging to the continuous spectrum: 
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In general, we can have an arbitrary number (Mr) of ra- 
diation modes, and we assume that they belong to the 
components with Mg + 1 < m < Alg + Mr- Note that 
solution (^, which is valid in the limit r„i — » 0, reduces 
to a simple plane-wave profile given by Eq. in the 
absence of bright components. When the radiation wave 
amplitudes are not small, both the radiation modes 
and the bright soliton profiles defined by Eq. (^) are dis- 
torted according to the nonlinear superposition principle. 
In order to find MSG solutions on a finite background, 
we extend the approach introduced in Ref. [Q. First, 
we multiply Eqs. (|^) by u*/{kj — ifx), add the complex 
conjugate, and sum over the fundamental soliton num- 
bers 1 < i < Mg. Comparing the resulting expression 
with Eq. (0), we obtain the following relation: 
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where the subscript m indicates a radiation mode. The 
second step is to define the combined rotational-scaling 
transformation which changes the nonlinear eigenmode 
amplitudes while preserving the total intensity profile up 
to a constant background. 
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Here / is the full intensity of the re-scaled solution 
Um{x,z), and lb = J2jl^M^+i\'''j\'^ ^^'^ background 
intensity. Then, from Eq. (ph, it immediately follows 
that the rescaling coefficients for the bright components 
should be defined as 
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while for the radiation modes we put Um = 1. Quite 
remarkably, this procedure can be used for both signs 
of nonlinearity (s = ±1). The propagation constants 
should also be modified accordingly, so that the resulting 
functions satisfy the original Eqs. (0). We finally obtain: 



(11) 



At this point, the derivation of the analytical solutions 
for the MSCs existing on top of several radiation modes 
is cornplete, and the component profiles are defined by 
Eqs. <^ and (|l]) together with Eqs. (|) and (0). All 
such solutions correspond to multi-parameter families 
which can be generated with the help of the rotation 
transformation (||). 

Let us now extend the analytical results to the case 
where the background is composed of a continuum set of 
radiation modes, i.e. Mr — > -|-oo. This case corresponds, 
for example, to spatial optical solitons excited by an in- 
coherent light source [^,0. Then, the plane waves in 
the background can be characterized by an angular dis- 
tribution function, -R(/i) (> 0), so that R{fi)dfj. is the 



wave intensity corresponding to the tangents of inclina- 
tion angles in the interval (/i,/i -I- dfi). Therefore, the 
full background intensity is lb ~ J-^ R{n)dii, and the 
scaling coefficients for bright components are 
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We note that, for a finite number of radiation modes, 
the distribution function can be written as R{fJ.) = 



fj,j), and then expression ( |12D reduces 



to Eq. (0). 

The above results are valid for both self-focusing (s = 
+1) and self-defocusing (s = —1) media. As follows from 
Eq. (|9|), the qualitative difference is that, in the former 
case we have bright complexes on a constant background 
while in the latter case dark dips are formed. 



IV. INTERACTION PROPERTIES OF MSCs 

The general analytical solutions can be used to obtain 
some important characteristics in a simple form. In par- 
ticular, we can use the fact that the outcome of a soliton 
collision, in terms of total intensity, is not affected by 
the presence of the background. Therefore, according to 
earlier results the interactions of MSCs are char- 
acterized by translational shifts of fundamental solitons 
after the collision, and these are given by 
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Here the summation involves the fundamental solitons 
which feature in the collisions with the soliton number 
j. When the colliding soliton number m comes from the 
right (i.e. has a larger x coordinate before the impact), 
we put fjm = +1 while we set fjm = — 1 when it comes 
from the left. The shifts differ for each soliton in an MSG. 
As a result, the intensity profile of an MSC changes after 
collisions. 



V. MODULATION OF BACKGROUND 
COMPONENTS 

According to the general relation (||), the intensity 
profile is uniquely determined by the eigenvalues of the 
bright fundamental solitons and the background intensity 
lb, but docs not depend on the the angular distribution 
of radiation waves. However, the total intensity of the 
soliton components, 
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and the intensity of the radiation modes, 
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both depend on the scahng coefficients Um defined in 
Eq. dl^). As follows from Eq. (p^, each fundamental 
soliton creates a hole in the background, and the corre- 
sponding modulation depth is proportional to the bright- 
dark coupling coefficient, given by the value (|J7mP — s). 
Quite interestingly, the radiation mode profiles are the 
same in self-focusing and self-defocusing media, provided 
the distribution function and soliton eigenvalues remain 
unchanged. 

However, there are some key differences between soli- 
tons in self-focusing and self-defocusing media. In the 
former case, a modulation of the background is compen- 
sated by the bright components having larger amplitudes 
(since \Um\'^ > s, and s = +1). On the other hand, in 
a self-defocusing medium, a dark soliton creates an ef- 
fective waveguide, which in turn can trap bright com- 
ponents. Such a self-trapping mechanism results in the 
limitation of the minimum dark soliton width. This hap- 
pens because the maximum intensity contrast is limited 
by the value of the background intensity. As a matter of 
fact, the limitation can be even stricter, since the max- 
imum modulation depth, ^A — maxxih — I)/Ib < 1, 
cannot always reach the value of 1. Then, according to 
Eq. (|^), the characteristic width corresponding to one 
fundamental soliton cannot exceed the value (Ailb)^^^^ ■ 
The actual limit is determined by solving the existence 
conditions, which follow from the requirement for the 
right-hand side of Eq. (^2|) to be non-negative, since, by 
definition, \Um\'^ > 0. It is interesting to note that these 
conditions involve only the individual wavenumbers of 
fundamental solitons, and they are automatically satis- 
fied for interacting solitons forming MSCs. 

Finally, we note that the radiation modes are charac- 
terized by a non-trivial phase modulation. For practical 
applications, it is especially important to know the phase 
jump, or the additional phase shift which appears due 
to the presence of bright fundamental solitons. Using 
Eqs. (^ and we find the following relation, 

■ I, 

ifi + k^' 

m=l ^ '"■ 

where (j) is the phase jump, and defines the inclination 
angle of the radiation mode. Then, the phase jump can 
be found as a sum over the phase shifts associated with 
individual fundamental solitons, 

/ \ 

'/'(m) = X! '^™('^) = X! 2arctan — I . (15) 

m=l m=l ^'-^ 

We see that the absolute values of the individual phase 
shifts are limited to tt. However, the total phase jump 
can become larger than tt if M > 1, i.e. if the MSG is 
composed of several fundamental solitons. 




FIG. 1. (a) Total intensity profiles for multisoliton 
complexes, and the corresponding decomposition between 
(b) bright and (c) dark components in a self-focusing medium. 
Each MSG consists of three fundamental solitons with 
Tm = 2,3,4; the MSG on the right is stationary (prn = 0), 
while the MSG on the left has a positive velocity (/im = 3). 
Angular width of background distribution is p = 1. 



VI. GENERAL RESULTS FOR A GAUSSIAN 
ANGULAR DISTRIBUTION 



Let us now analyse the features of the bright-dark de- 
composition for a Gaussian-type angular distribution, 

R{lJi)^h-^e~^'"/p\ (16) 

where p > is the characteristic angular width. Since 
the integral in Eq. (|l^) cannot be expressed in elemen- 
tary functions for arbitrary p, we first consider the lim- 
iting cases. Specifically, for a narrow angular distribu- 
tion, i.e. p <C 1, we have (|?7mP - s) ~ h/{rln + t^m), 
while for p ^ 1 (and p fim) we obtain (|?7mP — s) — 
Ib\/T^ / {prm) 0. Therefore, we expect that, for a fixed 
background intensity J;,, the modulation of the radiation 
waves should be reduced (i) for wider angular distribu- 
tions, i.e. larger p, and (ii) for MSGs having higher ve- 
locities \lJirn\- 

On the other hand, since the phase jump depends on 
the radiation mode wave number /i, the excitation of soli- 
tons can be more difficult in cases of wider angular spec- 
tra of radiation modes, i.e. larger p. Additionally, for a 
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moving MSG, i.e. when = const 7^ 0, the dependence 
4>{^i) becomes asymmetric unless R{p,—firn) — R{^^'m — ■ 




X 



FIG. 2. Intensity distributions for p = 7 in a self-focusing 
medium. Parameters and notation are the same as in Fig. hi. 



VII. BRIGHT MSG IN A SELF-FOGUSING 
MEDIUM 



that the lateral shift of the MSCs is relatively large. For 
example, it can easily be seen on the scale of Fig. ||. In 
contrast to single solitons, MSCs experience larger shifts 
in collisions, due to the multiple contributions from all 
the constituent fundamental solitons. 




FIG. 3. Interaction of two multisoliton complexes existing 
on a multi-component background in a self-focusing medium. 
Input profile corresponds to Fig. |l| 



For a self-focusing nonlinearity (s — +1), MSCs ex- 
ist in the form of bright localized waves having higher 
intensity than the background, as follows from Eq. (H). 
Examples of the total intensity profiles and the bright- 
dark mode decomposition are shown in Figs. ^ and g^for 
different values of p. The MSC on the right has zero ve- 
locity (with the corresponding = 0), while the other 
MSC (on the left) has a positive velocity; thus in the 
latter case, the corresponding background modulation is 
smaller, as predicted in Sec. VI. 

A collision between two MSCs is illustrated in Fig. ||. 
This example corresponds to the initial conditions shown 
in Fig. |l| A remarkable fact is that the total intensity 
profile does not depend on the value of p, provided that 
lb is preserved. The intensity profile for the collision will 
be the same for other values of p or for other distribu- 
tion functions. In these examples, the MSC actually has 
an intensity which is relatively small compared with the 
background level. 

Note that the shape of each MSC changes after a col- 
lision, for the reasons discussed in Sec. IV. In particular, 



a symmetric MSC becomes asymmetric after a collision 
(see also Ref. 0]). The presence of radiation does not 
influence this process. Another feature of a collision is 



VIII. DARK SOLITONS IN SELF-DEFOGUSING 
MEDIUM 



To describe MSCs on a background in a self-defocusing 
medium (s = —1), we first have to determine the ex- 
istence conditions, as outlined in Sec. 0. Considering 
the case of the Gaussian distribution given by Eq. (p^), 
we find that, for a narrow angular spectrum, in the 
lowest-order approximation, the existence condition is 
^"m + A'm — ^b- Therefore, the minimum soliton width. 



which is of order ; 



can be achieved if the soliton veloc- 



ity is zero. In the other limit where p 3> 1 and p 3> Pmi 
we have < /fj-y/Tr/Pi i-^- the minimum width increases 
linearly with an increase in p. Numerically-calculated 
existence regions are shown in Fig. ^ for two values of p 
in Eq. (p^). The figure clearly shows that the existence 
region is very similar to that in the case of a single compo- 
nent radiation field when p is relatively small. However, 
the existence regions become visibly different when p is 
large. 

The minimum soliton width versus p is shown in 
Fig. ||. This result shows that the distribution func- 
tion for the radiation field influences the properties of 
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an MSG, namely it changes the limiting parameters for 
the existence of an MSG, although the intensity profile 
of the MSG is not directly influenced by the properties 
of the radiation field. 



(a) 




-20 -10 10 



FIG. 4. Grey shading marks existence regions in the pa- 
rameter space of fundamental soliton eigenvalues (rm . Hm ) ■ 
The angular distribution function is given by Eq. (uM with 
lb ~ 40 and (a) p = 1 or (b) p = 7. Dashed lines correspond 
to the case of a single component background, when p — > 0. 
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FIG. 5. Dependence of the minimum characteristic soli- 
ton width on the parameter p for the angular distribution 
function of radiation waves given by Eq. ( |l6[ ) with = 40. 

Figure ^ shows an example of intensity distribution for 
the case of a self-defocusing medium (i.e. dark MSCs 
on a background). We have chosen the soliton eigen- 



values to be the same as in Figs. |^ and ||. According 
to our general expression the total intensity profiles 
in self-defocusing and self-focusing media are "mirror- 
images" relative to the level of the background. Even 
the radiation mode intensities coincide in these two cases, 
cf. Figs. |l](c) and ^(c). However, the bright component 
intensities arc different, as is clearly seen in Fig. |l|(b) 
and |^(b). This is a manifestation of the nontrivial na- 
ture of the nonlinear superposition of the solitons and 
the background components. 

Figure ^ shows the collision of two MSCs on a back- 
ground. Again, we can see that the nonstationary in- 
tensity profile created by the soliton interaction during 
collision is the "mirror image" of that for bright MSCs in 
a self-focusing medium, as shown in Fig. ^. The symme- 
try relation is mathematically exact. Correspondingly, 
the lateral shift is also governed by the same rules as 
those for a bright MSG. 

An important consequence is that the change of refrac- 
tive index induced by incoherent MSCs has exactly the 
same pattern in cases of self-focusing and self-defocusing 
media with Kerr-type nonlinearity. This remarkable fact 
can be used to design linear multimode waveguides and 
X-junctions with special properties. Moreover, the MSCs 
are controlled by many parameters (fundamental soliton 
eigenvalues), and therefore a wider variety of require- 
ments can be satisfied in comparison with X-junctions 
formed by a collision of two solitons. 




FIG. 6. Intensity distributions in a self-defocusinj 
medium. Parameters and notation are the same as in Fig. 11 
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FIG. 7. Interaction of two multisoliton complexes exist- 
ing on a multi-component background in a self-defocusing 
medium. Input profile corresponds to Fig. ^ 

IX. CONCLUSIONS 

In conclusion, we have obtained an exact solution 
for multisoliton complexes on top of a multi-component 
background composed of radiation waves in Kerr-type 
nonlinear media. We have identified similarities and dif- 
ferences between bright and dark MSCs which exist in 
self-focusing and self-defocusing media, respectively. In 
particular, we have found that the intensity profiles in 
these two cases are "mirror-images" relative to the level 
of the background, and that they depend only on the 
eigenvalues of the fundamental solitons. For example, 
the reshaping of MSCs after collisions is determined by 
the lateral shifts of the fundamental solitons and is not 
affected by background components. On the other hand, 
the width of dark solitons has a minimum, and this value 
depends strongly on the angular distribution of the ra- 
diation waves. We have derived analytical estimates of 
the key soliton characteristics for the case of a Gaussian 



angular distribution of radiation waves, and presented 
numerical examples illustrating the principal features of 
bright and dark MSCs. 
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